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We derive analytic merger rates for dark-matter haloes within the framework of the Ex- 
tended Press-Schechter (EPS) formalism. These rates become self-consistent within EPS once 
we realize that the typical merger in the limit of a small time-step involves more than two pro- 
genitors, contrary to the assumption of binary mergers adopted in earlier studies. We present 
a general method for computing merger rates that span the range of solutions permitted by 
the EPS conditional mass function, and focus on a specific solution that attempts to match 
the merger rates in iV-body simulations. The corrected EPS merger rates are more accurate 
than the earlier estimates of Lacey & Cole, by ~ 20% for major mergers and by up to a fac- 
tor of ~ 3 for minor mergers of mass ratio 1 : 10 4 . Based on the revised merger rates, we 
provide a new algorithm for constructing Monte-Carlo EPS merger trees, that could be useful 
in Semi-Analytic Modeling. We provide analytic expressions and plot numerical results for 
several quantities that are very useful in studies of galaxy formation. This includes (a) the 
rate of mergers of a given mass ratio per given final halo, (b) the fraction of mass added by 
mergers to a halo, and (c) the rate of mergers per given main progenitor. The creation and de- 
struction rates of haloes serve for a self-consistency check. Our method for computing merger 
rates can be applied to conditional mass functions beyond EPS, such as those obtained by the 
ellipsoidal collapse model or extracted from iV-body simulations. 

Key words: cosmology: theory — dark matter — galaxies: haloes — galaxies: formation — 
gravitation 
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1 INTRODUCTION 

The hierarchical clustering of dark matter is the key pro- 
cess in establishing the observed structure in the universe. 
Galaxies form inside the potential wells induced by the 
dark-matter distribution. The building blocks of this hier- 
archy are virialized collapsed gravitating systems in pres- 
sure equilibrium — the dark-matter haloes — characterized 
by their growth history, structure, and clustering. Although 
dark-matter dynamics is governed solely by the gravitational 
force, we are still far from a good quantitative understanding 
of its various features. 

The Press-Sche chter (PS) formalism 

dPress & Schechterl 1 19741) has been very useful in mod- 
eling the abundance of dark-matter haloes as a function of 
mass a nd ti me. It has been f urther developed bv lBond et al.l 
(fl99l and lLacev & Cold (fl99l hereafter LC93) to the 
Extended Press-Schechter (EPS) formalism, which provides 
at any time the mass function of progenitors of a halo of a 
given current mass. EPS has been a basic tool for under- 
standing the growth history of haloes, and it has been shown 
to grasp many of the key features of the buildup of haloes 
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in cosmological iV-body simulations (e.g. lLacev & Cold 
fl99ilCole etaill2008tlNeistein & Dekelll2008l) . While EPS 
has been used extensively for the last two decades, it still 
involves central open issues. One is the construction of 
self-consistent Monte-Carlo merger trees for Semi- Analytic 
Models of galaxy formation. The other is how to compute 
halo merger rates that will be consistent with the EPS mass 
function. 

While drawing the basic lines of the EPS theory, LC93 
worked out a formula for the merger rates of haloes. This 
formula has been popular in many applications, although 
it involves a problem. LC93 themselves noticed that their 
merger-rate formula has a problematic intrinsic asymmetry 
between progenitors of mas s M and Mp — M (whe re Mq is 
the descendant halo mass). ISheth & Pitman! ( 1 19971) realized 
that the LC93 assumption of binary mergers is not accurate 
when the power spectr um differs from a w hite noise (their 
discussion near eq. 27). lBenson et al.l (120051) interpreted this 
as an intrinsic inconsistency within the EPS formalism. We 
show below that the typical mergers have multiple progen- 
itors, more than two, even in the limit of a small time step. 
Adopting this correct limit, we obtain accurate EPS merger 
rates, which improve the LC93 estimates and are fully con- 
sistent with the EPS conditional mass function. The error 
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in the LC93 formula makes a significant difference for the 
number of merger events and for the fraction of halo mass 
added by mergers. 

Random realizations of merger trees that follow the 
EPS conditional mass function are widely used as the 
back-bone of semi-analytic modeling of galaxy forma- 
tion. Several different metho ds fo r constructing such 

trees h ave been proposed JCo le 1991; Kauffma nn & White] 
19931; ISheth & Lemsonl 119991; ISomerville & Kolattl 119991; 



Cole et al. 120001: Hiotelis & Popolo l 20061) . In most cases, 
these algorithms fail to recover th e EPS mass func- 
tion. I t seems that the algorithm of Kau ffmann & Whitel 
(119931) is the on l y one that is fully consistent with EPS. 
Shet h & Lemsonl d 19991) described an alternative that is also 
accurate, but it has not been developed into a detailed so- 
lution. One can indeed show that the EPS formalism per- 
mits many different types of merger trees that recover the 
EPS progenitor mass function. We provide below a new al- 
gorithm for constructing EPS merger trees based on our for- 
mula for merger rates. This algorithm does reproduce the 
EPS progenitor mass function, and it is chosen among the 
different solutions to be a good match to the merger trees 
extracted from cosmological A-body simulations. 

Empirical algorithms for generating merger trees that 
resemble the trees in cosmological A-body simulations have 
been p roposed by iParkinson et al.l d2008l) ; iNeistein & Dekell 
(12008b . These merger trees are for most parts better approx- 
imations to the TV-body results than any EPS-based tree. 
However, a correct EPS model has several useful benefits. 
For example, it allows very high mass resolution at low 
cost, it can be easily applied within any desired cosmologi- 
cal model, and it is self-consistent with the Press-Schechter 
halo abundance. On the other hand, the empirical algorithms 
mentioned above should always be verified and re-calibrated 
when used in a different cosmology or when applied at a 
different resolution. Analytic models in the spirit of EPS 
can serve us in understanding several open issues concern- 
ing the way haloes are identifi ed in A-body simulation s. For 
example, it has been noticed (Nei stein & De kel 2008, here- 
after ND08) that some of the non-Markov features in A- 
body merger trees may arise from the way haloes are de- 
fined. Indeed, the halo definition has become an open issue 
with the finding that the range of virial equilibrium in small 
haloes can extend well beyond the tra ditional "virial ra- 
dius" that is based o n spherical collapse dCuesta et al.ll2007l: 
iLudlow et al.ll2.Q08b . As part of this work we provide a gen- 
eral method for generating merger trees that follow any given 
conditional mass function. This mass function could be ei- 
ther based on sph erical collapse (i.e., EPS ), or arise from el- 
lipsoidal collapse dSheth & Tormenl2002l) . or extracted from 
A-body simulations. 

This paper is organized as follows. In $2] we present 
nomenclature, describe the limit of small time-steps, and 
prove the theorem concerning multiple progenitors. In £|3] 
we address different solutions for the EPS halo merger rates, 
and choose the solution that fits well the A-body results. In 
§4]we work out useful results for merger rates from our EPS 
formalism, and present them in practical formulae and in fig- 
ures. In ||5]we address the creation and destruction rates of 
haloes. In fj6]we describe a Monte-Carlo algorithm for con- 
structing EPS merger trees based on our adopted solution. In 
i|7]we summarize our results and discuss them. 



2 GENERAL ANALYSIS 



2.1 Definitions: PS and EPS 



In the EPS formalism, the natural dimensionless time vari- 
able is w(z) = 5 c (z)/ D(z), where D{z) is the cosmolog- 
ical linear growth rate of density fluctuations as a function 
of redshift z and 5 C ~ 1.69. The natural mass variable is 
S(M) = a 2 (Ad), the variance of the initial density fluctua- 
tion field, linearly extrapolated to z = 0, and smoothed using 
a window function that corresponds to a mass M. The reader 
is referred to ND08 for our specific way for computing these 
quantities. The cosmological model used here is defined by 
(Q A , O m , h, cr 8 ) = (0.75, 0.25, 0.73, 0.9), with the power 
spectrum specified in ND08. This model was adopted to en- 
able comparison with res ults extracted from th e Millennium 

cosmological simulation (ISpringel et alj|2005b. 

According to the EPS formalism dBond et al.l 119911 
LC93), the average number of progenitors in the mass in- 
terval [M, M + dM], which will merge into a descendant 
halo Mq after a time-step Aw, is given by 
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where AS = S(M) — S(Mo). We term the most massive 
progenitor in this time-step by Mi, the second most massive 
by M2, and so on. The probability that M is the mass of 
the i-th progenitor is termed Pj = Pi(Al\Mo, Aw). Conse- 
quently, the sum of all the P^'s equals dA/dA/: 



P tot (M|M ,Aw) 



dJV, 

— (M I M , Aw) 



Pi(M\M ,Au) 
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For brevity, we may sometimes omit the explicit dependence 
of Ptot and Pi on Mo and Aw. 

It is often useful to define a minimum halo mass, M m j n . 
Haloes with smaller masses are considered to be part of 
a smooth accretion component, encompassing a total mass 
M acc . 

We also need the total number density of haloes per unit 
mass per comoving volume, which is given by the Press- 
Schechter mass function: 
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where p is the present mean mass density of the universe. 

2.2 Number of Progenitors in a Small Time-Step 

Throughout this work, we appeal to the limit of a small time- 
step, Aw — > 0, relevant for the derivative with respect to 
"time", d/dw. For given Mo and M m - ln , the limit of a small 
time-step is defined here as Aw <C S(M — M min ) — S(M ). 
In this limit, and when M ^ Mo 
Ptot can be written as 



M m ; n , the probability 



P to t (M I M , Aw ^0) 



1 M Aw 
W M~(AS) 3 / 2 
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dM 



(4) 



Merger Rates of Dark-Matter Haloes 3 



10 



M =10 a 




10" 



M 



10 

. /M 



10" 



o 



Figure 1. The average number of progenitors given that the main progeni- 
tor mass is less than M max = Mo — M m ; n , as a function of M m j n /Mo. 
The three different curves are for different values of Mo as indicated 
(with units of h^ 1 Mq). Each of the curves is plotted only for M m ; n > 
10 6 (i _1 M0. The computation is done in the limit of a small time-step, 
Auj — » 0. Evidently, if the minimum mass is less than ~ 10~ 3 Mo, the 
number of progenitors is larger than two. This implies that the concept of 
binary mergers is highly inaccurate for low values of M m ; n /Mo. This con- 
clusion is valid independently of the value of Mo . 



after the exponent in eq. [TJ has been set to unity. Conse- 
quently, the "time" derivative of Ptot is simply 



dPtot(M|A/ ) 



1 M Q 



2^ M {AS) 3 / 2 



dS 



dM 



(5) 



We occasionally write d/dcj when it should formally be 
d/dAw, as both derivatives are the sam^B The above equa- 
tions are valid only for M ^ Mo — M m i n ; otherwise AS 
may also become infinitely small, such that ( Alu) 2 / AS does 
not vanish, and the exponent in eq. ([TJ) does not converge to 
unity. 

We now prove the theorem of multiple progenitors, 
claiming that according to EPS, the typical merger involves 
multiple progenitors rather than a binary merger even in the 
limit of a small time-step. Theorem: Given the EPS pro- 
genitor mass function of eq. Q}, with the CDM power- 
spectrum, in the range M min -C M and in the limit 
Auj — > 0, the average number of progenitors per merger 
event is greater than two. 

We first notice that the constraint of mass conservation, 
that the total mass in progenitors cannot exceed Mq, implies 
that events with M\ > M max , where M max = Mq — M m i n , 
cannot have any other progenitor with Mj > M m j n . There- 
fore, merger events between two or more progenitors above 
M m i n are limited to the cases where Mi < M max . 

Let TV be the number of progenitors with mass in the 
range [M min , M max ]. We first show that (N\Mi < M raax ) > 
2. If P(Mi < M) is the probability that Mi < M, then 
(N) = P(Mi < M max ) x (7V|Mi < M max ), because the 



1 We assume that Auj 
a fixed luq . 



luq and the derivative d /duj is computed at 



contribution of the other events is zero progenitor^ We thus 
obtain 



{N\M\<M n 



(N) 



P(M x <M max ) 

1- f™° P 1 (M)dM 
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When we calculate the integral in the denominator, we note 
that Pi can be replaced by P to t near Mo (see the discussion 
preceding eq. [9] below). As Alu — > 0, eq. implies that 
(N) vanishes in proportion to Aw, but P(Mi < A/ max ) also 
vanished, making the ratio converge to a finite value. 

Figure [TJ shows the average number of progenitors 
given that the main-progenitor mass is smaller than M max , 
(N\M\ < Mmax). This is in the limit of a small time-step 
and for different values of Mq. We see that this average is 
greater than two for any A/ m j n < 10~ 3 A/o. It increases with 
decreasing M m ; n to a value of ~ 10 for M m i n = 10~ 6 Mo. 
This proves that (N\Mi < M max ) > 2. 

Since each of the events with Mi < M max that are 
not mergers contributes to the conditional average of TV a 
value «C 1, the merger events, which are a subset of the 
Mi < M max events, must have on average even more pro- 
genitors than computed in eq. (0 and shown in Fig. [TJ We 
conclude that the assumption of binary mergers is invalid in 
EPS, even for Alu — > 0, once M m i n < 10 _3 Mo. This proves 
the theorem. 

If Mo is not that much larger than M m ; n , the range 
M min > 10 _2 A/ in Fig.Q] we obtain (iV|Mi < M max ) < 2. 
This implies that the average mass of the two progenitors 
does not sum up to Mo, namely the accretion component 
M acc contains a non-negligible fraction of the mass. 

Since the theorem of multiple progenitors has interest- 
ing implications on the formation of structure, it would be 
worthwhile to consider analytically the average number of 
progenitors in the idealized case where the power spectrum 
is a pure power law, S oc Af~ Q . Solving for (JV|Mi < 
M max ), we find that it is bigger than 2 for any < a < 1 
(once M m i n is small enough). For a = 1, the case of Poisson 
white noise, one can show that (iV|Mi < M max ) — > 2 when 
M m i n /Mo — » 0, in agreeme nt w ith the coagulation appr oach 
discussed bv lEpsteini ( 1 1 983b and lSheth & Pitman! d 19971) . For 
a > 1, the average number of progenitors never exceeds 2. 
We learn that the average number of progenitors per merger 
event depends on the shape of the power spectrum. In par- 
ticular, for power spectra that are relevant on galactic scales, 
a < 1, the average number of progenitors per merger are 
more than two. 

We note that the existence of multiple mergers in 



Strictly speaking, we should write (N) = P(M m i n < Mi < M max ) X 
(iV|M min < Mi < M max ). Using the fact that P(M min < Mi < 
M max ) ^ P(Mi < M max ) it is evident that all the results proved here 
using P(Mi < M max ) are lower limits on the accurate (TV |M m ; n < Mi < 
M max ). 

3 When computing the denominator one should use eq.fT] and not the ap- 
proximation of eq. f4] In the case where Af m i n <C Mo we can approx- 
imate Mo/M ~ 1 and the denominator is just Erf [Au;/\/2AS m ] ~ 
V / 2A"Au;/v / AS' m , where AS m = 5(M max ) - S(M ) and Auj -* 0. 
When Af m ; n is not small enough, the integral in the denominator can be 



computed only numerically. 
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the limit of small time-steps is already me ntioned in 
ISheth & Pitman! (Il997l) . lsheth & Lemsonl (Il999h added that 
for a general power spectrum, one can group the progeni- 
tors into sub-groups that merge like the progenitors of the 
Poisson-power-spectrum case. 

2.3 Merger rates 

One way to define a merger rate is as the probability for the 
i-th most massive progenitor to merge into the main pro- 
genitor within a time-step Alu. This is the joint probability 
for the two progenitor masses M\ and Mi, which we denote 
Pi j i(Mi, Mi\M , Alu). Note that P and P 3 can both have 
non-vanishing values at the same mass, so the probability 
for any progenitor with mass M s to merge with Mi is the 
sum 

P ltS {M u M s \M , Alo) = ^ P M (Mi, M S \M , Alu) . (7) 

i 

We learned in ^2.21 that there are typically several pro- 
genitors in each merger event even in the limit of a small 
time-step. To complicate matters even further, we note that 
P\.i+i is not necessarily smaller than P\.i- Still, for the pur- 
pose of estimating merger rates, we wish to approximate 
this multi-progenitor merging process as an instantaneous 
sequence of binary mergers. There is clearly no unique way 
to do that. We adopt here the assumption that each of the 
secondary progenitors (Mi,i > 1) merges with a halo of 
mass Mi, and ignore mergers among the secondary haloes 
themselves. The validity of this assumption can be tested in 
TV-body simulation^. This assumption makes sense when 
Mi is much more massive than the other progenitors. How- 
ever, in a case where Mi ~ M2 ^ M3, one might con- 
sider Ms merging with a halo of mass Mi + M% instead. 
We assume that this uncertainty in interpreting the multiple 
merger events does not translate to a significant error in our 
estimated average merger rate, but the actual estimate of this 
uncertainty is beyond the scope of the present paper. 

For any progenitor Mj, we define 
Pj| 1 (M I |Mi, Mo, Aw) to be the conditional probabil- 
ity to have Mi given that the main-progenitor mass is Mi. 
Then 

P 1A (M 1 ,M i \M ,Au)= (8) 
P^MilMi, M , Alu) • Pi(Mi|M , Alu) . 

Our approach for finding a solution for Pi.; starts with a so- 
lution for Pi, followed by a solution for P^. This is because 
Pi is determined robustly by EPS, with only a small, control- 
lable uncertainty over a limited mass range. 

The shape of Pi, the small freedom in it within EPS, and 
its effect on the av erage mass history of the main progenitor 
has been studied in lNeistein et alj d2006l) and can be summa- 
rized as follows: In the range Mi ^ Mo/2, Pi is identical to 
the known Ptot, because any progenitor in this mass range 
is by definition the main progenitor. For Mi slightly below 

4 In !0 we discuss the possible relation between a multiple merger event in 
EPS and a correlated sequence of binary mergers in an ./V-body simulation. 
If the time between mergers in the TV-body sequence is shorter than the time 
it takes the remnant halo to settle into its new potential, than our assumption 
might be reasonable. 
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Figure 2. Two possible solutions for the probability distribution of the main 
progenitor, Pi . The thick blue curve corresponds to Pt ot , properly normal- 
ized as indicated with Mo = 10 12 h^M® and Au) = 10~ 6 . The so- 
lutions for Pi differ only in the small tail at M < Mo/2. The shaded area 
marks the range over which the integral of Ptot equals unity; it ends at 
xi = M/Mq ~ 0.44. This is also our default definition of Pi, termed 
sharp tail. The dashed curve marks another possible tail, also correspond- 
ing to an integral of unity, which is linear in M, and thus termed linear tail. 
Note that we plot Ptot x Mo / Au as this curve is the same for all small 
Au;, in accord with eq.f4] 



Mo/ 2, there is a "tail" of non-vanishing probability, which 
could go to zero in many different ways. This is illustrated 
in Fig. |2] which shows two of the many possible solutions 
for this tail. Our default option is with a 'sharp tail', 

f P tot if Mi > xi M 
Pi(Mi|M ,Aco) = <^ . (9) 

I if Mi «C x x M 

The value of x\ is set by the requirement that the integral 
over Pi equals unity. It is ~ 0.44 for the cosmology and 
for the halo masses used her^El We note that the average 
mass history of the main-progenitor using this Pi can b e 
computed by the analytical formula of lNeistein et al.l d2006l) . 
Fig.|2]also shows an alternative solution where the Pi tail is 
linear in M. The freedom in the tail of Pi corresponds to 
an uncertainty of less than 8% in the average relative growth 
rate of the main progenitor, (dMi/dw)/Mi. 

Assuming a specific solution for Pi, the constraints for 
having a correct P^i are as follows: 

P(M|M ) = J P t |i(M 4 |Mi,M ) Pi(Mi|M ) dM 1 (10) 
P tot (M) = ^p(M) (11) 

i 

P(M X ,M 2 ,...) = if ^M;>M ( 12 ) 

i 



5 Using u = Iogio(Mrj) — 12 where Mo is in units of h Mq we can 
approximate xi = 7.118X 10 -5 « 3 +6.225x 10 _4 m 2 +0.0035«+0.444 
with an accuracy that is better than 0.05%. 
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The last condition is assuring mass conservation, where the 
total mass of all progenitors cannot be larger than Mq. 

For certain purposes, it will be helpful to define P s h as 
the sum over all P^. The constraint for P s |i is simply 

P tot (M) - Pi(M) = J P^iMlM^P^M^dM, . (13) 

Here mass conservation cannot be formulated as an explicit 
condition on P s |i because it does not contain information 
concerning the mutual distribution of multiple progenitors. 

While Pi is robustly determined in EPS, there is a 
great deal of freedom in P\^- This is because P^j is a two- 
dim ensional function wi th only one-dimensional constraints 
(e.g. lBenson et al.ll2005l) . We emphasize that this is true also 
for small time-steps. Hence there are many solutions for the 
desired EPS merger rates. In the next section we show sev- 
eral valid solutions of this sort. 



3 SPECIFIC SOLUTIONS 

Here we bring a general formalism for obtaining solutions 
Pi i and demonstrate the level of freedom allowed while 
obeying the EPS conditional mass functions. Given the ro- 
bust expression in eq. for P to t in the limit of a small time- 
step, the solutions presented below are valid for any value of 
Au> once it is small enough. 




Figure 3. Three solutions for Pi, 2, with a unique M.% for each M\. The 
solutions are derived here for M = 10 13 /i _1 Mq, Au; = 10 -6 ; they 
are practically the same for any smaller Aoj. The solid (blue) and dashed 
(red) curves are computed for the same Pi (the default sharp tail), and they 
differ only in the initial conditions (termed solutions / and II). The dotted- 
dashed (green) curve is obtained using the linear tail for Pi (solution ///). 
Note that the dashed and dotted-dashed lines have disconnected segments 
near (Mi , M 2 ) ~ (M /2, 0). The solid line is our default solution (I). A 
summary of these solutions can be found in tableff] 



3.1 Determining a unique set of Mi 's for a given Mi 

Our general solution is motivated by the merger-rate concept 
introduced by LC93. Assume that for any Mi we can choose 
a unique set of smaller progenitors {Mi}, so that each P^ 
is a delta function: 

P^iM^MuMo, Aco) = 6 [Mi - fi(Mi\M , Alu)} . (14) 

Here /j(Mi|Mo, Aim) associates a value of Mi to any Mi. 
We often write /j(Mi) where Mq and Aw are obvious from 
the context. Substituting P^ from eq. ( [Pfl ) in the constraint 
of eq. ( [Tol l, and integrating over Mi, we obtain a differential 
equation for /i(Mi): 

d/i(Mi) _ P(Mi) 
dMi Pi [/((Mi)] ' ^ 

where / is assumed to be a monotonically decreasing func- 
tion of Mi. Thus, the solution for fi(M\) is determined by 
Pi, Pi, and a certain initial condition M o = /i(Mi o). This 
differential equation is to be integrated numerically to obtain 
a solution for /j(Mi). Note, in contrast, that LC93 adopted 
the inaccurate assumption (Mi) = Mo — Mi, failing to 
allow for the additional progenitors beyond M2. 

We start, for example, with P2J1, using our default 
sharp-tail solution for Pi as in eq. (|9). Given this Pi, we 
try to set P2 = Ptot — Pi, which simply equals P to t in the 
range M < xiMo. A solution for /2(Mi) can now be ob- 
tained for a given initial condition. Our first choice, which 
we term "Solution I", is 

(Mi,o, M 2 , ) = (x 1 M , Xl M ) . (16) 

This ensures that M2 approaches Mi as the latter ob- 
tains its minimum value xiMq. Solution / is shown in 



Fig. [3] We also plot the solution for the initial condition 
(M^o, M2.0) = (Mo — X1M0, xi Mo), termed solution 77. 
As is evident from the figure, although both solutions have 
the same Pi, they have quite different values of Pi, 2- Fig- 
ure [3] also shows Solution 7/7, which is based on Pi with 
the "linear tail" shown in Fig. [2] A summary of these three 
solutions is listed in tableQ] It should be noted that solutions 
II and 777 include a small range of Mi values near Mo/2 
that is not connected to M2 through M2 = /2(Mi). For Mi 
in this range we cannot use any value of A/2 that was previ- 
ously associated with Mi through f2- We can either choose 
Ma = 0, or treat it similarly to P3, as will be explained be- 
low. 

Luckily, the condition for mass conservation is almost 
fully obeyed by each of the three solutions for fi{M{) 
above. This is implied by the fact that the curves for Pi, 2 
seem to always lie below the line Mi + M 2 = Mq. How- 
ever, a closer look shows that this constraint is violated for 
Mi > 0.99M . In this range f 2 (M l ) > M - Mi for 
all the solutions presented here, and we cannot adopt the 
M2 that solves the differential equation. Instead, we en- 
force M2 = Mo — Mi, which makes the distribution P2 
differ slightly from P to t- This result is expected based on 
the multiple-progenitor theorem of i)2.2| requiring more than 
two progenitors for reproducing Ptot- 

Fig. |4] shows Ptot and P2 for small M values. The ef- 
fect discussed above leads to P2 < P to t at M2 < O.OlMo, 
meaning that additional progenitors are needed in order to 
obtain an accurate fit for Ptot- 

Next we should address Pj for for the i-th progenitor, 
i > 2. In what follows we use as an example the solution / 
of Pi, 2, and the procedure can be easily generalized to deal 
with the other solutions. For i > 2 we define 
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Figure 4. The probability distribution for small progenitors according to 
our algorithm, eqs. {T7} and {T8), using M = 10 13 h~ 1 M Q and Aui = 
10 — 6 . The solid line is Ptot . The dashed (blue) curve is P2. It deviates from 
Ptot for M < 6 X 10 — 3 X Mo - The Dashed-dotted curves correspond to 
Pi for i > 2. They equal Ptot — Pi- The vertical bars mark the limits for 
each Pi, termed A/high,! and Mi OWj i. 



Table 1. The characteristics of the EPS solutions for P\,i discussed in i]3.1l 
The threes solutions assume the delta-function form for P211, eq. 1141 . 



Solution 


Pi tail 


Pi, 2 Initial conditions 




Sharp 


(x 1 M ,x 1 M ) 


11 


Sharp 


(xiMo, 1 - 11M0) 


III 


Linear 


(M /2, Mq/2) 



f P tot (M) - P 2 (M) Mi OWii <M<M higV 
P(M) = <^ . (17) 

t otherwise 

where A/high, ?: = A/i OWii _i for i > 3 and A/high, 3 is the max- 
imum A/ for which P2 < Ptot - The value of M\ ow ^ is set by 
the condition of mass conservation: each solution fi(Mi) is 
defined up to the point where Mq = A/i + ^ fi{Mi). The 
initial condition is thus 

(Mi, , M <l0 ) = (aJiMo, A^ higM ) , (18) 

and the set of Pi we obtain is given in fig.|4j 

To summarize, our solution for the merger rate is 

P 1 ,i(M 1 ,M i \M ,Au>)= (19) 
Pi(Mi|M 0j Alu) S [Ah - /i(Mi|M , Aw)] . 

In the limit of a small time-step, /; does not depend on Alu 
and Pi is given by eq. (0]). 

3.2 Comparison with TV-body results 

We now wish to compare the merger-rates from our EPS 
analysis to merger-rates that were extracted from TV -body 
simulations. We remind the reader that when using TV -body 
simulations with small time-steps, the merger-rates suffer 
from inconsistencies due to non-Markov features (ND08), so 



any Markov model, as the one implied by EPS, will have de- 
viations at small time-steps. A more fair comparison should 
be done against the model derived by ND08, which should 
be close to the optimum Markov fit to the simulations. 

In Fig. [5] we show results of our solutions I and 
II against TV-bod y simulations (the Millennium run, 
ISpringel et al.l2005h . In order to test our solutions we use the 
merger-tree algorithm as described in section|6]below. Fig. 7 
of ND08 indicates that the EPS merger rates found here 
do resemble relatively well the merger rates of the Markov 
model that fits the simulation in ND08. At bigger time-steps, 
we see that although the general contour shape is similar, the 
average mass of the second progenitor is slightly smaller in 
E PS than in the simulatio n. This is also evident in the results 
of iParkinson et al.l d2008l) . who compared a different set of 
TV-body merger trees with a binary-merger model for EPS 
trees a la LC93. On the other hand, merger tr e es con structed 
using the algorithm of ISomerville & Kolattl d 19991) have a 
significantly lower mass for M2, as p ointed out in ND08. 
We fin d that the algorithm proposed bv lKauffmann & White! 
( 1 19931) produces EPS trees that match the TV-body trees at a 
level comparable to our EPS solution /, though it may not be 
as useful as our algorithm in generating a statistical sample 
of merger trees and in allowing analytic estimates. 

When comparing Fig. 7 of ND08 to Fig. here, it 
seems that the most important difference lies in the shape 
of the main-progenitor distribution. We recall that according 
to ND08, for trees extracted from TV-body simulations, this 
distribution is log-normal in S. On the other hand, in EPS 
this distribution is given by eq. (0, which has quite a differ- 
ent shape (see also Fig.|2]i. 

Figure|5]also compares solutions / and //, showing that 
solution / is somewhat closer to the TV-body results. At the 
smaller time-step, solution II has slightly higher values of 
Mi than the simulation, and it also has an isolated peak near 
(Mi, M2) = (Mo/2, 0), with no parallel trace in the sim- 
ulation. At Auj = 1.7, solution II shows bigger deviations 
in the masses of Mi and M2. Based on these findings, we 
adopt solution / as our default option for Pi, 2- However, one 
should bear in mind that each of the three solutions discussed 
above is an example of a solution that is fully consistent with 
the EPS conditional mass functions. 

3.3 A More Realistic Model? 

The solution in terms of delta-functions, eq. (fl4] i. is moti- 
vated by the work of LC93 and by results from TV-body sim- 
ulations. We find that the Pi 2 extracted from the Millen- 
nium simulation indeed approaches a narrow function when 
Aw — > 0. Nonetheless, it should be noted that the delta- 
function solution is not the only possible solution for EPS 
even when Aw — ► 0. We do find other EPS solutions with 
a broad P2|i- The delta-function treatment is simple, though 
it has its limitations. For the finite time-steps used, the ac- 
tual width of the distribution in the TV-body simulation is fi- 
nite, not zero. With the optimal time-step for reconstructing 
merger trees, Auj ~ 0.1 (ND08), the delta-function solution 
is accurate within EPS, but it is not such a good approxima- 
tion to the TV-body merger trees. 

A different approach might be to seek a solution that 
can be used with any time-step Alu in a self-consistent way, 
namely it should keep the same when using k time-steps of 
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Figure 5. The joint probability of the two most massive progenitors, Pi, 2, 
for haloes of mass 2 X 10 13 h~ Mq. The plots refer to two differ- 
ent time-steps, Aw = 0.1 and 1.7. The contour levels are at Pi 2 — 
5, 10,30 Mq 2 . Upper panel: The solid (red) contours refer to results of the 
EPS algorithm I, obtained from generated random realization with intrin- 
sic Ao; = 0.001. Lower panel: Same for solution II. The corresponding 
results from the Millennium TV-body simulation are shown for comparison 
as dashed (blue) contours (previously presented in ND08, Fig. 7). 



Awi or one time-step of Aoj = k x Au>i. Motivated by 
ND08, we try 

P 2 | 1 (M 2 |Mi, M , Alu) = P x {M 2 \M a ,Au) , (20) 

where M a = f 2 (M±) as defined in eq. (fl4l . This solu- 
tion is fully consistent with P to t for small enough Aw as 
it approaches a delta function. However, for big time-steps 
it shows some deviations from the theoretical Ptot> depend- 
ing on the specific solution adopted for Pi. 2- This solution 
is not practical for our applications because it does not fit 
accurately the shape of Pi j2 as obtained from many small 
time-steps of solution I. We mention it here because it is 
close to solution 77 even for big time-steps. 

For completeness, the explicit expression for the merger 
rate in this case is 



P lj2 (M 1 ,M 2 |M ,Aw) 



1 M M, 



(Ac) 2 



2tt MiM 2 (ASiASa) 1 - 5 



(21) 



100 



3 

■D 




M 10 



Figure 6. The number of merger events with mass ratio > r, per unit "time" 
du>, for a given final halo mass Mo (in units of h" 1 Mq). The value of r 
associated with each pair of curves is indicated. The solid curves describe 
the results of our EPS model. The dashed curves are the results of the LC93 
formula, assuming that the main progenitor is more massive than 0.5Afo 
- M - Mi. 



and tf c 



xexp 



(Aw) 2 / 1 



1 



2 \ASi AS 2 



dS(Mi) 



dMi 



dS(M 2 



dA/ 2 



where A Si =S(Mi) — S(Mq) and AS 2 = S(M 2 )-S(M a ) 



4 PRACTICAL APPLICATIONS 

We next implement the method outlined above for EPS 
merger rates, specifically solution I, to compute several 
quantities concerning the clustering of dark-matter haloes, 
which are of practical interest in the studies of galaxy for- 
mation. 



4.1 Major and Minor Merger Rates for a Given Halo M 

We first compute the probability that a halo of mass Mq has 
undergone within the last time-step Aui a merger event that 
includes the main-progenitor M\ and another progenitor of 
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mass Mi > rM\ (i ^ 2). Using our definition for merger 
rate, section 1231 this can be written as 



diV big (r,U ) 
duj 



(22) 



— E 



P hi {M 1 ,M i \M , Aw)dMidMi 



We use eq. ( fT9l l in order to integrate over Mj and obtain 

diV big (r,M ) / dPi(Mi|Af ' 



duj 



= E 



fi(M 1 )>M 1 r 



doj 



dMi . (23) 



Eq. <(5j provides the derivative of Pi with respect to uj for any 
r > 0. This rate is independent of redshift as it is expressed 
in terms of the self-invariant time variable uj. When needed 
in units of time, one should multiply the above expression 
by uj. A useful approximation for u> (from ND08) is 

u = -0.0470 [1 + z + 0.1(l + z) -1 - 25 ] 2 ' 5 h 73 Gyr" 1 ,(24) 

where /173 is the Hubble constant in units of 73 km s _1 . This 
approximation is valid for the ACDM cosmology used here, 
with (O ro ,0 A ) = (0.25,0.75), to better than 0.5% at all 
redshifts. 

Figure [6] shows results for diVbig(r, Mo)/dco. For ex- 
ample, we read that dA big (0.3, 10 12 h^M^/du ~ 0.65, 
which means that a halo of mass 10 12 ft. _1 M© has under- 
gone on average 0.65 major mergers of r > 0.3 per unit of 
uj. Multiplying by Co at z = gives 0.04 major mergers per 
Gyr. At z = 3 it yields ~ 1 such mergers. The number of 
minor mergers, with 10~ 4 < r < 0.3, is drastically higher; 
a 10 12 /i _1 Mq halo has ~ 10 such minor mergers per Gyr 
at z = 0, and ~ 250 such events per Gyr at z = 3. 

Figure|6]also shows the number of merger events as de- 
rived from the formula of LC93, and assuming that the main 
progenitor is more massive than Mo/2. The LC93 approach 
is interpreted here as /^ c (A/i) = Mq ~ M%. The error due 
to their assumption is ~ 20% for major mergers, and it be- 
comes as large as a factor of ~ 3 at r ~ 10~ 4 . We emphasize 
that this is true for our default solution /. It is possible that 
another EPS solution may be somewhat closer to the LC93 
results, but the discrepancy of the LC93 estimates for minor 
mergers is likely to remain large. 



4.2 Growth Rate of a Halo M due to Major Mergers 

As a second example we compute the average mass fraction 
added to a halo by merger events with mass ratio greater than 

r, 



dF hig (r,M ) 



du) 
d_ 

duj 



(25) 



Mi 



— y / Pi,i(Mi, Mi\Mo, Auj)^ L dM 1 dM i 



As before, we can simplify the expression to 

cLF big (r,M ) 



dw 

E 



(26) 



dPi(M 1 |M )/ i (Mi) 



fi(M 1 )>M 1 r 



du> 



Mn 



dMi 



3 
jo 



5> 




0.01 



M 10 



Figure 7. The mass fraction added to a halo of mass Mo by mergers with 
progenitors of mass ratio > r, The values of r are the same as in Fig. [6] The 
solid curves describe our EPS model. The dashed curves refer to LC93; they 
are plotted only for the extreme values of r. 



Results for dF b i g (r-, Mo)/dw are shown in Fig. [7] As an 
example, dP big (0.3, 10 12 hr l M Q )/duj reads - 0.2. This 
means that a halo of mass 10 12 h~ 1 M Q has gained on av- 
erage ~ 20% of its mass by major mergers per unit of uj. 
Multiplying by uj we get a growth rate of 1% per Gyr by ma- 
jor mergers at z = 0, and ~ 30% at z = 3. For this quantity 
the LC93 assumption leads to similar errors of < 20% for 
all mass ratios r. 



4.3 Merger Rates for a given Mi 

The number of haloes of mass M s that will merge with a 
halo of a given mass Mi (M s < Mi) within the time-step 
Aoj, with no restriction on the descendant mass Mq, is 



dQ(M s \Mi,z) 
doj 



(27) 



Ml/X1 dPi, s {Mi, M S \M , Aoj) <t>{M ,z) 
duj <(>(M 1 ,z) 



dM 



M ± +M s 
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where (j)(M, z) is the Press-Schechter average comoving 
number density of haloes of mass M at redshift z, eq. (|3j. 
Using Pn from eq. dT9b we get 



&Q(M 8 \Mx,z) 



(28) 



E 



dPi(Mi|A/ o , t )0(M o ,^) 

&LQ <f>(Ml,z) 



dfi 



dMn 



where Mo,i are the values of Mq for which M s = 
fi(Mi\M 0ii ). Note that the derivative d/,/dM is with re- 
spect to Mo rather than Mi . 

Figure [8] shows results for dQ(M s \Mi, z)/duj. Unlike 
the other quantities discussed above, Q does depend explic- 
itly on redshift z, through the dependence of the sum in 
eq. (EHT l on z. Nevertheless, for major mergers (high r) this 
sum consists of only one term, so the z dependence can be 
scaled out. Note also that Q is not a smooth function, due to 
the fact that fi(Mi) are always defined for an Mi value that 
is smaller than some threshold, in order to conserve mass 
(see the discussion after eq. (fT8l). Figure|8]displays in com- 
parison the results of the LC93 formula, showing deviations 
of ~ 30% for major mergers, which become as large as a 
factor of ~ 3 for a small mass ratio of r ~ 10~ 4 . 



5 CREATION AND DESTRUCTION RATES OF HALOES 

In this section, we address the merger processes through the 
creation and destruction rates of haloes. We show that the 
results obtained here are consistent with the Press-Schechter 
mass function, as they should be. This can also serve as a 
sani ty check for validatin g the values of Q computed above. 

iBenson et al.l (12005b have attempted to use the Smolu- 
chowski coagulation equation for computing halo merger 
rates. This equation evaluates the change in the number den- 
sity of haloes due to the competing processes of halo cre- 
ation and destruction as a result of mergers. Their formulae 
assume that each halo is formed by a binary merger event, 
so halo formation is modeled as a 2-progenitor process. Ac- 
cording to the multiple progenitor theorem proved above, 
this approach cannot yield correct results. A formulation of 
the Smoluchowski coagulation equation should be therefore 
replaced with a different equation that takes into account 
multiple mergers a nd accretion mass. Co nsequently, the dis- 
crepancy found by IBenson et al. I (120051) in the merger rate 
formula of LC93 is not a discrepancy in the EPS formalism - 
it simply reflects the inaccuracy introduced in the LC93 for- 
mula by the assumption of binary merg ers. Despite this built- 
in error, the numerical merger rates by IBenson et al.l (12005b 
seem to be consistent with the Press-Schechter mass func- 
tion. This could be the result of using a mass grid cell size 
of Mo/179, which scales with mass, while we found that 
multiple mergers become relevant only for M < 10~ 3 Afo. 

The EPS formalism is, by construction, fully consis- 
tent with the Press-Schechter mass function. This can be ex- 
pressed by 



<j){M,u>- 



-Aw) = \ 

JM 



^(Mo, W )Ptot(M|M ,Aco)dMo .(29) 



This equation indicates that any merger rate that is consis- 
tent with P to t must be consistent with the way varies in 
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Figure 8. The number of mergers of mass M s = rM\ with a given halo 
of mass Mi, in an infinitesimal time-step duj. The results are plotted for 
z = 0. The solid curves are results of the EPS solution / given in [[3] The 
dashed curves follow the formula of LC93. 



time. This implies that the merger rates that were evaluated 
here should predict the correct rate of change of when im- 
plemented using halo creation and destruction terms. 

Taking into account the multiple progenitors and the ac- 
creted mass, the time derivative of (j>(M, z) is connected to 
Q via the equation 



lim 



1 



d<f>(M, z) _ 

duj AAf- > oTAw-»o Alo 2AM 

M+AM i-M— AM 

0(M O , z)dM / Pi (Mi\M , Auj) dAfi 

M—AM JO 

pM+AM 

0(M O , z)dM / Pi (Mi | M , Auj) dMi 

M+AM 

dQ{M\Mi,z) 



M 



dw 



/ M—AM 

'-4>(M 1 ,z)dM 1 



(30) 



The first term corresponds to the creation of new haloes in- 
side the mass bin [M — AM, M + AM] as arising from 
the main-progenitor growth rate. The second term computes 
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the number of haloes that leave this bin for the same reason. 
We note that each of these terms diverges for small AM, but 
their sum remains constant. The third term that involves Q is 
the number of haloes that leave the mass bin by merging with 
bigger haloes. We have verified that this formula yields self- 
consistent results by computing it term by term. However, 
due to the numerical limitations of computing Q in only dis- 
crete points, we get an accuracy that is on the order of few 
percents in the integral of Q. 



6 A MONTE-CARLO ALGORITHM FOR EPS MERGER 
TREES 

Using the specific analytical solution obtained in $3] one can 
construct full merger trees. The algorithm is conceptually 
simple, and can be summarized as follows: 

• Define a reference halo with mass Mq at c^o- 

• Choose a time-step Auj (not necessarily small). 

• Draw a random main-progenitor mass Mi, using the 
distribution Pi (Mi |M , Alu). 

• Compute the value of Mj {i ^ 2) using Mj = 
fi(Mi\Mo, Aw), for every value of i, until the desired mass 
resolution is achieved. 

• Repeat the above procedure for each progenitor Mj, 
where Mo is replaced by Mj. 

This general algorithm can be used with any variant of the 
solutions presented in $3] An advantage of this algorithm 
is that all the tree quantities can be computed analytically. 
Another advantage over other algorithms is that its accuracy 
within EPS is in principle unlimited — it solely depends on 
the accuracy of the /j used. The algorithm can be applied 
with time-steps that are not small, but the procedure is sim- 
pler when using small time-steps so P t ot is linear in Auj. 

Figure [9] shows results from EPS merger-tree realiza- 
tions using our algorithm based on solution /. These results 
demonstrate the high accuracy of the generated trees. 



7 SUMMARY AND DISCUSSION 

We presented a rigorous method for computing dark-matter 
merger rates and merger trees that obey the halo progeni- 
tor mass function of the EPS formalism at an y redshift. This 
corre c ts apparent inconsis tencies within EPS dLacev & Cold 
119931: iBenson et al.ll2005l) . Our method conserves mass, in 
the sense that the sum of the progenitor masses does not ex- 
ceed the mass of the product halo. This method translates the 
problem of constructing merger trees to solving a differential 
equation. Different choices of initial conditions correspond 
to different types of merger trees. This method enabled us to 
span the set of solutions for merger rates within EPS, and to 
pick up a specific solution whose merger trees are a good fit 
to TV-body results. The same method can be implemented 
with any conditional mass function beyond EPS, e.g., as 
extracted from TV-body simulations or from an ellipsoidal- 
collapse model. 

Our main result is an accurate derivation for the merger 
rate of dark-matter haloes, w hich differs from the classical 
result o flLacev & Cold 0993). This is due to our finding that 
within the EPS formalism, a merger event typically involves 



many progenitors in a time-step, even when this time-step is 
infinitely small, as opposed to the binary mergers assumed 
in previous works. Our c o rrecte d results differ from those 
derived by lLacev & Cold d 19931) especially in the number 
of minor merger events, while other quantities deviate only 
at the level of 20%. We compute a few useful variants of 
the merger-rate formula, such as the number of mergers for 
a given descendant halo, the mass fraction added by merg- 
ers, and the merger rate per progenitor halo. These exam- 
ples span many applications for galaxy formation models. 
We also verified that the merger rates derived here are fully 
consistent with the evolution of the Press-Schechter mass 
function, in terms of counting the creation and destruction 
of haloes within the coagulation equation. 

We have shown that the merger rates derived here fit 
the resul ts of TV-body si mulat ions better than the early re- 
sults of Lacev & Cold (119931) . However, as discussed in 
iNeistein & Dekell (120081) . the merger rates from TV-body 
simulations may suffer from intrinsic inconsistencies at the 
level of a few tens of percents due to non-Markov ef- 
fects. Keeping this in mind, it is tempting to compare our 
EPS merger rates with other studie s of merger-rates ex- 
tracted from TV-bod y simulations (e.g. IFakhouri & Mal2007t 
IStewart et al.ll2.Q07b . For example, it is likely that our EPS 
results are in better agreement with the TV-body resul ts than 
the EPS results presented by iFakhouri & Mai d2007l) : their 
EPS merger rates are underestimates at low mass ratio of 
r ~ 10~ 3 . This better agreement is similar to what we fin d 
here based on the merger-rates of lNeistein & Dekel d2008l) . 

The concept of multiple mergers in the limit of small 
time-steps, proven here to be valid in EPS, deserves further 
attention. Recent studies indicate that this might be true in 
TV-body simulations when the time-steps used are finite (e.g. 
IFakhouri & Mall2007t INeistein & Dekelll2008l) . However, in 
an TV-body system, every merger event can be broken into 
a sequence of binary mergers once the time step is short 
enough. This implies that the TV-body system is not a pure 
Markov process — the binary mergers are a non-Markov 
feature. This non-Markov feature reflects a correlation be- 
tween the successive mergers. A correlation of this sort may 
be introduced, for example, by the progenitors being part 
of a cosmic-web filament feeding a bigger halo, where they 
merge in as a coherent group. When we impose a Markov 
model to describe the TV-body mergers, i.e. a model that 
ignores any correlations, this correlated sequence of binary 
mergers is forced to appear as a multiple merger. It would 
be interesting to verify this interpretation of the relation be- 
tween the EPS and TV-body mergers by testing whether most 
of the TV-body merger events are indeed part of a correlated 
sequence. 

If multiple-merger events are a common phenomena in 
EPS, then the merger rates are not defined in a unique way, 
as the counting method by which progenitors are ordered 
to merge with each other may affect the merger-rate results. 
Here we have chosen a simplified approach, where all the 
progenitors are assumed to merge with the most massive pro- 
genitor and not with one another. Clearly, other methods of 
counting may be applied. This issue may be examined in de- 
tail using an TV-body simulation, where the multiple mergers 
can be broken into a sequence of binary events once the time 
steps are made small enough. It should be noted that the con- 
ditional mass function of progenitors as extracted from TV- 
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body sim ulations cannot be recon structed by a Markov pro- 
cess (see iNeistein & Deke]||2008l) . This means that there is 
no accurate expression for this mass function at small time- 
steps that can reproduce the mass function at high redshift, 
na mely separated fr om the present by a large time-step (but 
see ICole et al.ll2008L for an approximation). Further effort is 
needed in order to understand this issue in TV-body simula- 
tions. 

A conditional mass function that is based on the ellip- 
soidal collapse model has been used re cently for generat- 
ing merger tr ees dMoreno & S heth 2007) and for computing 
merger rates ( tZhang et al.ll2008l) 7The use of the ellipsoidal 
model is partly motivated by its earlier success, over the 
spherical model used by Press-Schechter, in reproducing the 
(un-c onditional) mass functio n of haloes in TV-body simula- 
tions (Sh eth & Tormenl2002l) . The method developed in this 
paper can be easily generalized to utilize the ellipsoidal col- 
lapse model. The results should be compared to our EPS pre- 
dictions and to the TV-body results. As a first step, it should 
be interesting to evaluate the level of accuracy in previous 
studies due to the binary-merger assumption by computing 
the average number of progenitors per merger event. 

The algorithm we provide for generating merger trees 
has several advantages as follows: 

• It is fully consistent with the EPS conditional mass 
function of progenitors. 

• The relevant statistics can be described analytically, in- 
cluding those concerning the main-progenitor history and 
the merger rates. 

• This algorithm was chosen, out of the many options that 
are consistent with EPS, to provide best fit to TV-body simu- 
lations. 

• The constructed merger trees conserve mass, in the 
sense that the total mass in progenitors does not exceed the 
descendant halo mass. 

These are significant im provements over previous algo- 
rithms that follow EPS Cold Il99lt iKauffmann & White! 



1991 ISheth & Lemsonl [19991: ISomerville & Kolattl [19991: 



Cole et al EoooT Hiotelis & P0P0I0II2OO6IC which makes the 
new algorithm a useful tool for analytic and semi-analytic 
modeling of galaxy formation. Still, the non-EPS algorithms 
that are empirically t u ned to match TV-body simulations 
dParkinson et al] 120081: INeistein & Dekell 120081) may have 
advantages in certain cases where the accuracy is important. 
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Figure Al. The number of merger events with mass ratio > r, per unit 
"time" du, for a given final halo mass Afrj (in units of 1i~ 1 Mq j ). The 
values of M are 10 8 , 10 1( \ 10 12 , 10 14 h -1 Af s . The solid curves 
describe the results of our EPS model while the dashed curves are the results 
of the LC93 formula. 



Figure A2. The mass fraction added to a halo of mass Mo by mergers 
with progenitors of mass ratio > r, The values of Mo are the same as in 
Fig. ED The solid curves describe our EPS model and the dashed curves 
refer to LC93. 
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APPENDIX A: MORE RESULTS 

This Appendix is a supplement to §4] for the benefit of prac- 
titioners who desire to read out numerical values for merger 
rates from the figures. Figures 6-8 of Sj4]present merger-rate 
quantities as a function of halo mass for different given val- 
ues of mass ratio r. Here we plot the same merger rate quan- 
tities as a function of r for different values of mass. This way 
of presenting the merger rates emphasizes their simple scal- 
ing with halo mass and highlights the trends at small values 
of r. 
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r 

Figure A3. The number of mergers of mass M B = rM\ with a given halo 
of mass Mi, in an infinitesimal time-step da;. The results are plotted for 
z = 0, and for Aft = 10 8 , 10 10 , 10 12 , 10 14 h^M®. The solid 
curves follow the EPS solution / given in ij3] and the dashed curves are 
obtained from the formula of LC93. 



